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Abstract. 

We study the interaction of dark-bright sohtons in two component Bose-Einstein 
condensates by suitably tailoring the trap potential, atomic scattering length and 
atom gain or loss. We show that the coupled Gross-Pitaevskii (GP) equation can 
be mapped onto the Manakov model. An interesting class of matter wave solitons and 
their interaction are identified with time independent and periodically modulated trap 
potentials, which can be experimentally realized in two component condensates. These 
include periodic collapse and revival of solitons, and snake-like matter wave solitons 
as well as different kinds of two soliton interactions. 
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The advent of solitons in ultra cold quantum gases has provided a new insight into 
the understanding of the localized wave packets that travel over long distances without 
attenuation. In this context, recent explorations of solitons in Bose- Einstein condensates 
(BECs) have paved the way for new developments in manipulating coherent matter 
waves for applications, including atom interferometry, coherent atom transport and 
quantum information processing or quantum computation. In particular, experimental 
observations of matter wave solitons of the dark [Ill2], and bright [3llll[5] types in BECs 
have attracted a great deal of attention in connection with the dynamics of nonlinear 
matter waves, including soliton propagation P [7], vortex dynamics [8], interference 
patterns [9] and domain walls in binary BECs [TO] . 

In this connection, recently much effort has been given to the study of matter 
wave solitons in BECs with time varying control parameters such as (i) variation of 
atomic scattering length which can be achieved through Feshbach resonance p!Tl [T2| 
[T3l [HI [15], (ii) inclusion of appropriate time dependent gain or loss term which can 
be phenomenologically incorporated to account for the interaction of atomic cloud or 
thermal cloud, and (iii) periodic modulation of trap frequencies ^6J. Bright solitons 
created in the experiments are themselves condensates and propagated over much larger 
distances than dark solitons which, on the other hand, can only exist as notches or holes 
within the condensates [Il[2|- Recently, dark solitons, their oscillations and interaction 
have been demonstrated in experiments with single component BECs [TTl [T8l [19] . Many 
studies on the bright/dark soliton formation and propagation of attractive/repulsive 
BECs have been focussed mainly on the single species systems. 

Multi-component generalization of the soliton dynamics is very natural in the 
context of atomic BECs because of the several ways to create such systems, for example 
as mixtures with two different atomic species/hyperfine states [20|,|2T1|22|, and as internal 
degrees of freedom liberated under an optical trap and atom-molecule BECs [23]. The 
multicomponent BECs, far from being a trivial extension of the single component 
ones, present novel and fundamentally different scenarios for their ground states and 
excitations [2H |25|. Matter wave solitons in multicomponent BECs hold promise for a 
number of applications, including the multi-channel signals and their switching, coherent 
storage and processing of optical fields. 

Since multi-component BECs are of greater interest for the aforesaid reasons, in this 
letter we investigate the interaction of dark and bright solitons in two component BECs. 
In the context of cold atomic gases, the two vector components which evolve under the 
Gross-Pitaevskii equation are the macroscopic wave functions of Bose condensed atoms 
in two different internal states, which we shall denote as |1) and |2). If one considers 
condensates of, for instance, ^^Na and ^''Rb atoms, the nonlinear interactions are due 
to elastic s-wave scattering among the atoms, and are effectively repulsive (positive 
scattering length) for both the systems in which multi-component condensates have 
been realized [26]. Very recently, there has been a tremendous interest in studying 
the dynamics of two-component Bose-Einstein condensates coupled to the environment 
using both experimental and theoretical means [271 ISHl [291 SOI (SH [32] . Here we bring 
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out exact bright-dark solitons in repulsively interacting two-component BECs coupled 
with the thermal clouds, which introduce gain/loss of condensate atoms. In particular, 
by considering both time independent and periodically modulated trap potentials, we 
show fascinating interactions of matter wave solitons. 

For the above purpose, we consider a trapped BEC with two components, where the 
dynamics takes place only in one dimension due to the strong trap confinement in the 
transverse direction. We have also included the interaction of external thermal clouds 
which is described by gain or loss term. The properties of BECs that are prepared in two 
hyperfine states can be described at sufficiently low temperatures by the dimensionless 
form of two coupled CP equations as [33] 



^j. J = 1,2. (1) 



k=l 

2 

Here V{x,t) = lQ'^{t)x^ is the external potential, fi^(t) = ^ , Ux is the trap frequency 
in the axial direction, uj± is the radial trap frequency {fl'^{t) < for expulsive potential 
and > for confining potential), n is the chemical potential, R{t) = as{t) is 

the s-wave scattering length, is the Bohr radius and gij = 1 and 7(t) = T{t) is 
the gain/loss term, which is the phenomenologically incorporated interaction of thermal 
cloud [Hll [35l |36] . Here we have considered the case of a two-component condensate 
with ^''Rb atoms prepared in two different hyperfine states for which the atomic masses, 
intra- and inter-component atomic scattering lengths are equal as described in Ref . [17] . 
Note that in equation ([T]), the variable x actually represents — , where a± 



and similarly t stands for u±t. 

In equation ([T]) when the gain/loss term j{t) is positive, it leads to the mechanism 
of loading external atoms (thermal clouds) into the BEC by optical pumping while 
j(t) < describes a BEC that is continuously depleted (loss) of atoms. If the condensate 
is fed by a surrounding thermal cloud, then the condensate undergoes an appropriate 
growth/collapse. In actual BEC experiments, for example, with ^Li and ^^Rb atoms, 
the atomic scattering length can be varied by suitably tuning the external magnetic field 
through the Feshbach resonance as [3], HI [15] 



where a° is scattering length of condensed atoms, B{t) is the external time varying 
magnetic field, Bq is the resonance magnetic field and A is the resonance width. Similar 
tuning of scattering length should be possible for two component condensates as well, 
for example in the case of hyperfine states of ^^Rb. 

Using the following point transformation [371 EHl IHHl HQ] 

^jix,t)=A,{x,t)q,iX,T), J = 1,2, (3) 

equation ([1]) can be reduced to the set of two coupled nonhnear Schrodinger equations 
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(2CNLS) of the form 



.dqj 1 d qj 
'df ^ 



2. 2 



^^ = -^gjli + (lj2^9M, J = 1,2, (4) 

k=l 



where Aj{x,t) = ro^2R{t) exp[i{9{x,t) + ^ijt) + / 7(t)dt], R{t) = R{t) exp [J -f{t)dt], 
9 = -hx^ + 2cirlRx - 2c\rl J R^dt, X = V2roRx - 2^201 / Ffdt, T = 2r^ J R^dt 
(ci, tq: constants), and R{t) and VL{t) have to satisfy the condition 

which is a Riccati type equation for Equation (jlj) is the so called defocussing 
Manakov system, which exhibits interesting one, two and A^-soliton solutions of bright- 
bright, dark-dark and dark-bright types [HI |42l |43] . Concentrating for the present only 
on dark-bright solitons, from the solutions of equation (jl]), one can straightforwardly 
construct the one, two and dark-bright soliton solutions for equation ([1]), provided Rt, 
7(t) and VL{t) satisfy equation We may note here that because of the complicated 
transformation involved above, the resultant soliton parameters are no longer constants 
but are functions of x and t. Therefore, the resultant soliton solutions may be considered 
as generalized solitons and not simple standard solitons. 

One-soliton dynamics:- The dark-bright components of the one-soliton solution of the 
defocussing Manakov system (jlj) can be given as 

where 77 = kX + i{^K^ — t'^)T, k = a -\- ib, p = k — ic, x = + 1^*)'^ — 1 j • The 
parameters a and b correspond to the amplitude and velocity of the soliton envelope, 
respectively, and c refers to the phase. 

Now using the transformation ([3]), the corresponding dark-bright one-soliton 
solution of the two coupled GP equation can be written as 

ijj{x,t) = \J 2R{t) qj{X, T) e^(e(-'*)+M*)+/ '^'^^ j = 1,2. (8) 

Two soliton dynamics:- The dark-bright two-soliton solution of the defocussing 
Manakov system (jlj) can be written as 



q,{X,T) = Le-^-^lcy^+r^]Tf^ _ y2 ^x^.e^^-^^*^ + P-lElj e^^^^i^^^^^'^ ) (9) 
d V Pk P1P2 



\j,k=l jtk=l 



<l2{X,T) = -\}^ e"^ - > ^ u,,xi,X2j e^^^^^^^' | , (10) 
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where 

2 

j,fc=i 

Vj = i^jX + i{\K] - t'^)T, Kj = aj + ibj, pj = Kj - ic, Xjk = + '^fc)^ ~ l) 

'^ifc = (/^i — f^kY (— — I" 1 )) / = XiiX22ki2Xi2p- The dark-bright two sohton solution 
of the two coupled GP equation ([1]) can be again represented by equations ([8]) with qi 
and q2 of the forms (jH]) and ffTUl) . The procedure can be extended to A^-soliton solutions 
also. However, we do not present their forms here. Depending on the form of the trap 
potential, gain/loss and interatomic interaction novel type of dark-bright matter wave 
solitons can be deduced using the above forms (l8l)- (fT0ll . In the following, we demonstrate 
them for two simple trap potentials. For the other choices, results will be presented 
elsewhere. In the present study we fixed the trap parameters similar to that used in a 
recent experiment on dark-bright solitons in two component ^''Rb condensates [17j . 

Time independent trap potential:- First let us consider, as an example, the case 
of time independent expulsive parabolic trap potential, Q'^{t) = —fig for which 
the integrablity condition ([5]) gives R{t) = sech(fio^ + The amplitude of the 




Figure 1. (Color online) Choice of the atomic scattering length as{t) and gain or 
loss term r{t) as a function of time that exhibits dark-bright solitons for (a) time 
independent expulsive trap potential, V{x) — —^Hqx'^ and (b) Periodic modulated 

. , , ,[ 4 + 10wcos(wi + (5) + 3t^^[l + cos^(a;< + (5)1 1 , 

trap potential, V(xA) = < 1 ^ — } . 

\ 4[l + wcos(wt + (5)]' J 



wave packet is given by ^ 2R{t)e^ '^^^'^ 1'^'^^ . Different types of soliton solutions can 
be constructed for suitably chosen gain term 7(t). In figure [U^a), we plot the gain 
r(t) = a;j_7(t) = 27ra;_Lf2o tanh(fio^ + '^)) and the corresponding choice of atomic 
scattering length as{t) = \aBR{t) = ^as sech.'^{ilot -\- 6), which may be realized by 
tuning the external magnetic field as 

aO A 
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We point out here that such a form of scattering length has been reahzed in ^Li and ®^Rb 
atoms [H m |15]. Figure [2]( a) shows the dark-bright one sohton solution for the above 
gain term where the amplitude of the wave packet remains constant. On the other hand, 
if we choose a periodic gain term, 7(t) = f2' sin(f2't), the amplitude shows oscillatory 
behaviour leading to periodic collapse and revival phenomena, see figure [2|^b). 




Figure 2. Dynamics of dark (left) and bright (right) components of the one soliton 
in an expulsive trap potential for (a) 7 — fio tanh(r2ot) and (b) 7 = 51' sin(r2't). The 
parameters are fixed at a = 0.7, h = 0.5, ro = l/v^, ci = c = 0, 17' = 2.5 and 
r^o = 0.001. 



We have also investigated different types of two soliton interactions for the gain 
term ^{t) = r2otanh[fio^] by suitably choosing the parameters ai, 61, 02, &2 and Qq in 
equations (Q and ffTOl) . Figure [3t^a) shows the interaction of dark and bright two solitons 
with strong beating effect without interaction while figure ^h) illustrates interaction 
of two solitons with beating effect. Figure [3](c) shows the soliton interactions without 
beating. These type of soliton interactions are well known in the fiber optics context 



Periodically modulated potential:- If we choose R(t) = 1 + a; cos{ujt + S) and 7(t) = 

— P " ^ . — ^-r-T, a; < 1 we get R(t) = \/l -\- u cosiut + S) and the integrability 

2[1 + uj cos{ujt + o)\ 

A-^- /m. • n2u\ 2/1 4: + 10ujcos{ut + 5) + 3u^[l + cos^{ujt + 6)]\ , . , 

condition @ gives Q (t) = < 1 ^ ^ — ^ ^ > which 

is a periodically modulated trap. For the above case, we sketch the gain T(t) = 

uj_i_uj'^ sinicut -\- S) , , 1. , . r . . 1 1 / \ 

— P 7 rrr and the corresponding choice of atomic scattering length asitj = 

2[1 + ucos{ujt + d)\ 

^ [1 + a;cos(a;t + 6)] in figure [U^b), which may be reahzed by periodically tuning the 
external magnetic field as 
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Figure 3. Contour plots of two soliton interactions sliowing dark (left) and bright 
(right) components in the presence of expulsive trap potential, (a) Strong beating of 
two soliton bound states without interaction for ai = 0.7, bi = 0, a2 = 0.77, 62 = 0, 
r^o = 0.001, (b) soliton interaction with beating for oi — 0.7, hi = —0.01, 02 = 0.9, 
&2 = 0.03, rif) — 0.01, and (c) elastic collision of two solitons for ai = 0.7, 61 = 0.3, 
02 — 0.7, &2 — —0.3, flo — 0.05. Here x axis represents time in nis with t £ (0, 200) 
for (a), t e (0,300) for (b), and t G (0,40) for (c). y axis corresponds to the axial 
coordinate in /im, with x £ (—6,10) for (a) and (b), and x E (—10,10) for (c). The 
other parameters are tq = 1/V2 and ci = c = 0. 



(a) (b) 




Figure 4. Dynamics of dark component (a) and bright component (b) of snake-like 
one soliton in the presence of periodic modulated harmonic potential. The parameters 
are a = 0.3, b = 0.2, ro = 1, ci = c = 0, (5 = and = 0.4. 



= ^0 + 1 hi ^ f ,^,,y (12) 
— ^ciB [l + co cos{cut + d)\ 

Figure 111(a) shows the snake-hke effect of the one sohton solution for periodically 

modulated trap potential with R(t) = \1 -\- u cos(ujt + 6)] and 7(f) = " ^ , — ^—rr. 

Similar snake effect has been demonstrated recently in scalar coupled nonautonomous 
NLS equations in BEC [52] and in the context of optical solitons jSj. Next we analyse 
different types of two soliton interactions for the periodically modulated potential and 
suitable choice of other parameters. Figure O^a) shows the dark and bright non- 
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interacting moving two solitons again with strong beating effect in bound state, wfiere 
tfie velocities are zero {pi = 62 = 0) for ai = 0.7, 02 = 0.8. Figure [5](b) shows non- 
interacting moving two solitons again with strong beating effect in unbound state for 
tti = 0.7, a2 = 0.8, 61 = 62 = 0.1. Figure [5](c) depicts the interacting two solitons with 
beating effect for ai = 0.7, 02 = 0.85, bi = 0.1, 62 = 0.12. Finally figure [St^d) depicts the 
two soliton interaction without beating effect for ai = 0.7, 02 = 0.85, bi = 0.1, 62 = 0.17. 




Figure 5. Contour plots of matter wave soliton interactions of dark (left) and bright 
(right) components in a periodically modulated trap potential (a) noninteracting two 
soliton with strong beating effect in bound state for ai ~ 0.7, 02 = 0.8, 61 = 62 = 0, 
(b) noninteracting two soliton with strong beating effect in bound state for oi = 0.7, 
02 = 0.8, 5i = 62 = 0.1, (c) two soliton interactions with beating effect for ai = 0.7, 
02 = 0.85, bi ~ 0.1, 62 ~ 0.12, (d) two soliton interactions without beating effect for 
ai — 0.7, a2 = 0.85, bi — 0.1, 62 = 0.17. The other parameters are ro = 1, ci = c = 0, 
uj — 0.1 and S — —5.0. Here x axis represents time in ms with t G (0, 300) and y axis 
corresponds to the axial coordinate in /^m, with x G (—10, 30). 

Recent experimental studies show the above types of dark and dark-bright soliton 
oscillations and interactions in single and two component BECs [T71 [T8| [T9]. However 
these experiments do not consider the interaction of thermal cloud or time dependent 
interatomic interaction. Our above studies clearly suggest the possibility of observing 
them in experiments. 

In summary, we have investigated the exact dark-bright one and two soliton 
solutions of the two-component BECs with time varying parameters such as s-wave 
scattering length and gain/loss term. On mapping the two coupled CP equation onto 
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the coupled NLS equation under certain conditions, we have deduced different kinds of 
dark-bright one sohton sohitions and interaction of two solitons for time independent and 
periodically modulated trap potentials. The present study provides an understanding 
of the possible mechanism for soliton excitations in multi-component BECs. These 
excitations can be realized in experiments by suitable control of time dependent trap 
parameters, atomic interaction and interaction with thermal cloud. 
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